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Abstract. The dynamical shear instability is an important mixing process in 
the advanced stages of the evolution of massive stars. We calculated different 
models of 15 Mq with an initial rotational velocity Ui n i = 300 km/s to in- 
vestigate its efficiency We found that the dynamical shear instability has a 
timescale shorter than Oxygen burning timescale and that it slightly enlarges 
the convective zones and smoothens the omega gradients throughout the evolu- 
tion. However, its effect is too localized to slow down the core of the star. 



1. Theory &: computer model 

The stability criterion for dynamical shear instability is the Richardson criterion: 
Ri = ( d Jjg z y2 > \ = where U is the horizontal velocity, z the vertical coor- 
dinate and iV 2 the Brunt- Vaisala frequency. Using this criterion, the timescale 
of the instability is the dynamical timescale and the instability is called dynam- 
ical shear instability (Endal & Sofia 1978). If heat losses are accounted for in 
Ri (Maeder 1997), the timescale is the thermal timescale and in that case the 
instability is called secular shear instability. The critical value, Ri c = 1/4, is 
used by most authors as the limit for the occurrence of the dynamical shear. 
However, recent studies (Canuto 2002; Briiggen & Hillebrandt 2001) show that 
turbulence may occur as long as Ri < ~ 1. 

Different formulae for the corresponding diffusion coefficient, D, are used 
at the present time (Zahn 1992; Maeder 1997; Heger et al 2000; Briiggen & 
Hillebrandt 2001). The following dynamical shear diffusion coefficient suggested 
by J. -P. Zahn is used in this study: 



where r is the mean radius of the zone where the instability occurs, AQ is the 
variation of O over the zone and Ar is the extent of the zone. The zone is the 
reunion of consecutive shells where Ri < Ri c . This is valid if the Peclet number, 



The computer model used is the Geneva evolutionary code. Modifications 
(to be described in a future paper) have been made to study the advanced stages 
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of the evolution of massive stars. The models are calculated for 3 stars of 15 
M Q at Z Q with v m i = 300 km/s. Convective stability is determined by the 
Schwarzschild criterion. The overshooting parameter, d ovcr /Hp =0.1 for El- 
and He-burning and afterwards. In the first model, the dynamical shear was 
not included. Then one model was calculated with Ri c = 1/4 and one with 
Ri c = 1. The calculations reached shell O-burning. Note that the computer 
model includes secular shear and meridional circulation for any calculation. 

2. Results &: discussion 

The characteristic timescale of the dynamical shear is very short (a fraction of a 
year) throughout evolution when using equation (1). We obtain diffusion coeffi- 
cients between 10 12 and 10 14 cm 2 /s. This is usually 1 or 2 orders of magnitude 
larger than the dynamical shear diffusion coefficients of Briiggen & Hillebrandt 
(2001) or Heger et al (2000). 

However, the extent of the unstable zone is very small, a few 10~ 3 M 
so the shear mainly smoothens the sharp il-gradients but does not transport 
angular momentum over long distances. The total angular momentum of our 
models, at the beginning and at the end of the calculations are of the same order 
of magnitude as the values obtained by Heger et al (2000) in their models E15 
and E15B. 

Although dynamical shear strengthens convection, especially the He-burning 
convective shell, the structure and the convective zones are similar between 
the model without dynamical shear and the one with dynamical shear and 
Ri c = 1/4. Concerning the Richardson criterion, we see that there is a dif- 
ference during Ne-burning between the model using Ri c = 1/4 and Ri c = 1. 
Indeed, using the latter, there is another large C-burning convective shell while 
the central convective core is smaller. One could argue that our diffusion coeffi- 
cient is too high but we obtained a similar trend using the diffusion coefficient 
of Briiggen & Hillebrandt (2001) with Ri c = 1. 

The well known discrepancy between the observed and predicted pulsar 
rotation periods might be reduced with the interaction between rotation and 
magnetic field (see contributions by Heger and Spruit in this volume) or with 
a formula linking dynamical and secular shear in the advanced stages of the 
evolution of massive stars. 
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